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Measurement of the cosmic-ray antiproton spectrum at solar minimum with a 
long-duration balloon flight over Antarctica 
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The energy spectrum of cosmic-ray antiprotons (p's) from 0.17 to 3.5 GeV has been measured using 
7886 p's detected by BESS-Polar II during a long-duration flight over Antarctica near solar minimum 
in December 2007 and January 2008. This shows good consistency with secondary p calculations. 
Cosmologically primary p's have been investigated by comparing measured and calculated p spectra. 
BESS-Polar II data show no evidence of primary p's from evaporation of primordial black holes. 



Precise measurement of the cosmic-ray antiproton (p) 
spectrum is crucial to investigations of conditions in the 
early universe and cosmic-ray propagation. Most cosmic- 
ray p's are produced by interactions of cosmic-ray nuclei 
vifith the interstellar gas. The energy spectrum of these 
"secondary" p's peaks near 2 GeV, decreasing sharply be- 
low and above due to the kinematics of p production and 
to the local interstellar (LIS) proton spectrum. The sec- 
ondary p's offer a unique probe |l|-t3j of cosmic-ray prop- 
agation and solar modulation. Cosmologically "primary" 
sources have also been suggested, including the annihi- 
lation of dark-matter particles and the evaporation of 
primordial black holes (PBH) by Havsfking radiation (J). 

Small PBHs, formed in the early Universe by initial 
density fluctuations, phase transitions, or the collapse of 
cosmic strings, might have a significant evaporation rate 
at the current age of the Universe and could_ contribute 
to the measured p spectrum at lovif energies [5[ . Because 
the predicted LIS PBH p spectrum peaks at ~ 150 MeV, 
this would be strongly influenced by solar modulation, so 
a search is most sensitive at solar minimum |6| . 

BESS95-I-97 showed that the p spectrum peaks around 
2 GeV [7|, and measurements by BESS and other exper- 
iments have shown that p's are predominantly secondary 
[8]. However, the low-energy p spectrum measured by 
BESS95+97 at the previous solar minimum was slightly 
flatter than predicted by secondary models. Although 
this suggested the possible presence of primary p's, the 
large statistical error of the BESS95-I-97 data did not 
allow a firm conclusion. BESS-Polar |8l4l3| was devel- 
oped to evaluate the possibility of excess low-energy p 
flux, with unprecedented precision, using long-duration 
solar-minimum flights over Antarctica. BESS-Polar I 
flew in December 2004 [H-li^ and BESS-Polar II [H] 



flew near solar minimum in December 2007 and January 
2008. Here, we report measurements of cosmic-ray p's 
from 0.17 GeV to 3.5 GeV by BESS-Polar II and dis- 
cuss the implications for secondary models and possible 
primary sources. 

BESS-Polar is a high-resolution magnetic-rigidity spec- 
trometer. A uniform fleld of 0.8 T is produced in a thin 
superconducting solenoid fllled with drift-chamber track- 
ing detectors. Particle trajectories are determined by flt- 
ting up to 52 hit points with a resolution of ^ 140/im in 
the bending plane, giving a magnetic-rigidity (= Pc/Ze) 
resolution of 0.4% at 1 GV and an overall maximum 
detectable rigidity (MDR) of 240 GV. Upper (UTOF) 
and lower (LTOF) scintillator hodoscopes measure time- 
of-flight (TOE) and dE/dx and provide the event trig- 
ger. For BESS-Polar II p measurements, the acceptance 
is 0.23 m^sr. TOE resolution between the UTOF and 
LTOF is 120 ps, giving a /3~'^ resolution of 2.5%. A 
threshold- type Cherenkov counter (ACC), using a silica 
aerogel radiator with optical index n = 1.03, rejects e~ 
and /Li~ backgrounds by a factor of 6100 to identify p's up 
to 3.5 GeV [ii]. A thin scintiflator middle-TOF (MTOF) 
on the lower surface of the solenoid bore detects low- 
energy particles that cannot penetrate the magnet wall. 
TOE resolution between the UTOF and MTOF is 320 
ps. In the present analysis, the MTOF was used to inde- 
pendently verify the procedure for eliminating interacting 
upward-going protons that could mimic low-energy p's. 

BESS-Polar II was launched on December 23, 2007, 
from Williams Field, near the US McMurdo Station in 
Antarctica, observing for 24.5 days with the magnet en- 
ergized. The float altitude was 34 km to 38 km (residual 
air of 5.8 g/cm^ on average), and the cutoff rigidity was 
below 0.5 GV. 4.7 x 10^ events were acquired with no 



TABLE I. p flux at the top 
observed p's and estimated 



of atmosphere with statistical (first) and systematic (second) errors, 
background events. The mean energy for each range was calculated 
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inflight event selection as 13.6 terabytes of data. 

In flight, most detectors and instrument systems op- 
erated well, with expected performance. Although the 
central tracker exhibited high-voltage fluctuations, nor- 
mal resolution was preserved for more than 90% of the 
observation time by using algorithms that calibrate the 
tracker over short time intervals and depend on its high- 
voltage state. Two TOF PMTs with high- voltage control 
problems were turned off, one on a UTOF paddle (of 10) 
and one on an LTOF paddle (of 12). Requiring two good 
PMTs on each paddle reduced acceptance ~ 20%. 

Analysis was performed as described in Ref. 
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FIG. 1. The /3~^ versus rigidity plot and p mass selection 
band after dE/dx and ACC cuts. For clarity, only 1 in 300 
positive-rigidity events is shown so only a few e+ or n^ can 
be seen. The lowest energy p's are shown in the inset flgure. 



same selection criteria were applied for p's and protons 
because they behave similarly in the symmetric configu- 
ration of BESS-Polar, except for deflection direction. 

Figure [T] shows f3~^ versus rigidity plots for events sur- 
viving dE/dx and ACC cuts. A clean, narrow band of 
7886 p's mirrors the protons. The calculated e~ and fi~ 
background is 0.0%, 1.0%, and 2.3% in the 0.2-1.0 GeV, 
1.0-2.0 GeV, and 2.0-3.5 GeV energy bands. Other back- 
grounds, such as albedo, mismeasured positive-rigidity 
particles, and re-entrant albedo, were negligible. 

The differential flux of p's at the top of atmosphere 
($toa) integrated over dE can be expressed as: 

$TOAd^ - (A^TOI - N,tmos)/e^n-/{Sn ■ Tn™) (1) 

A^TOI - {Np - NBG)/{edet ■ £non-int) (2) 

where Tuve is live time, and Np and Nbg are numbers 
of observed p candidates and expected background par- 
ticles. For the present analysis T^vo — 1286460 seconds. 
The effective geometric acceptance, including noninter- 
action efficiency {Sfl • enon-int), was calculated using 
GEANT3 as 0.133 ± 0.011 m^sr at 0.2 GeV and 0.159 ± 
0.008 m^sr at 2.0 GeV, with errors estimated from differ- 
ences relative to GEANT4. The detection efficiency for 
p's (£dct) was calculated using a noninteracting proton 
sample as 81.4 ± 0.1 % at 0.2 GeV and 60.0 ± 0.2 % at 
2.0 GeV. To obtain $toa, corrections were applied for p 
survival probability [l9| in the residual atmosphere (sair) 
and estimated atmospheric p production (iVatmos)- Eair 
was estimated as 85-6 ± 2.0 % at 0.2 GeV and 89.8 ± 
2.0 % at 2.0 GeV. iVatmos, 17.6 ± 3.2 % of the detected 
p's at 0.2 GeV and 27.6 ± 5.0 % at 2.0 GeV, was cal- 
culated by solving simultaneous transport equations [19| 
with adjusted interaction length (A) and tertiary produc- 
tion [20|. The uncertainty in this calculation is 18.1% 
(= (5.0%2(air depth) -K8.9%2(A)-hl5.0%2(tertiary))i/2). 
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FIG. 2. Solar minimum BESS-Polar II, BESS95+97 and 
PAMELA TOA p fluxes and secondary model calculations.. 

Table U gives the flux of p's at TOA from 0.17 to 3.5 
GeV with the statistical (first) and systematic (second) 
errors. The dominant systematics are atmospheric sub- 
traction and detection efficiency. A rapid change in ef- 
ficiency due to the ACC veto increases the systematic 
uncertainty in the two highest bins. 

Fig. [2] shows the BESS-Polar II p spectrum with 
BESS95-H97 and PAMELA [2l[ measurements and solar- 
minimum secondary calculations [y, |22h26| . Curve 1 uses 
Mitsui et al. [6|, \2^ data with force-field modulation 
of 600 MV from the best fit to the BESS-Polar II pro- 
ton spectrum. Curve 2 was generated by interpolating 
model calculations supplied by Bieber et al. 23] for neg- 
ative solar magnetic field polarity (A<0). The tilt angle 
of 15°(A<0) is the best fit to the BESS-Polar II pro- 
ton data. Curve 6 is the published A>0 solar- minimum 
calculation [23| for comparison to the BESS95-I-97 mea- 



surements. Curves 3 [2J| and 4 [25| are also published 



solar-minimum calculations. Curve 5 was generated us- 
ing the GALPROP model ^ with 600 MV force-field 
modulation. Improved statistical precision of the mea- 
sured p flux results from 14 and 30 times more events 
below 1 GeV than BESS95+97 and PAMELA, respec- 
tively. The BESS-Polar II and PAMELA spectra gen- 
erally agree in shape, but differ in absolute flux. The 
weighted mean difference, with combined uncertainties, is 
14 ± 5%, calculated near 2 GeV to reduce modulation ef- 
fects. Both are consistent with solar-minimum secondary 
calculations. Neither exhibits the flattening at low ener- 
gies found by BESS95-I-97, although the differences are 
statistically small. 

The evident differences among the calculations shown 
in Fig. [2] arise from several factors that can affect the 
normalization or shape of the spectrum: (1) definition of 
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FIG. 3. Comparison of p flux shapes with secondary calcu- 
lations normalized to BESS-Polar II flux at 2 GeV. Sensitiv- 
ity to uncertainty in the force-fleld modulation parameter is 
shown by the lower shaded band. The small sensitivity of 
secondary p drift calculations at solar minimum to tilt angle 
is illustrated by the upper shaded band. 



the primary proton and helium spectra, (2) incomplete 
knowledge of nuclear physics in propagation, (3) param- 
eters and models of propagation in the Galaxy, and (4) 
modulation in the heliosphere. Variation in the absolute 
fluxes of interstellar protons and helium, for instance, af- 
fects the absolute flux of p's, but not the spectral shape. 

Precise measurement of the low-energy p spectrum by 
BESS-Polar II allows secondary flux calculations to be 
evaluated by comparing observed and predicted spectral 
shapes, as shown in Fig. [31 The calculations are normal- 
ized to BESS-Polar II at 2 GeV to focus on their shapes. 
The calculated spectra and data points are also multi- 
plied by E'^^ to emphasize differences at low energies. 
The observed data are not normalized. Chi-square (x^) 
calculated with BESS-Polar II data and the normalized 
secondary p calculations in Fig. |3]are 0.61 (1), 0.61 (2), 
1.32 (3), 1.70 (4), 0.67 (5). The shape variation from 
uncertainty in the level of solar modulation is illustrated 
by the lower shaded band, calculated with the Mitsui et 
al. model [y, |2^ and modulation parameters of 500 MV 
(X^=0.81) and 700 MV (x^=0.52). The smah sensitivity 
of drift calculations to tilt angle is shown by the upper 
shaded band using A<0 Bieber et al. data [23| at 10° and 
20° (interpolated). In both cases, the change in spectral 
shape is small compared to differences arising from prop- 
agation models, because of the peaked shape of the LIS 
secondary p spectrum. BESS-Polar II results are more 
consistent with models (curves 1, 2, and 5) without low- 
energy p's from tertiary interactions (curve 3) or a soft 
spectrum from diffusive reacceleration (curve 4). 

The likelihood of primary p's from PBH evaporation 
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FIG. 4. (Top) Possible primary p fluxes from PBH evapora- 
tion calculated for BESS-Polar II (A) and BESS95-I-97 (B) 
by fitting differences of the measured spectra from the Mitsui 
secondary p spectrum. (Bottom) PBH evaporation rate (TV) 
distributions. Values of 7?. < are non-physical. 



can be quantified by a model-dependent evaporation rate 
(7^) determined by fitting a PBH model spectrum to the 
difference of a secondary calculation from the measured 
flux. TZ is positive (physical) only if the measured flux 
exceeds the secondary prediction. To avoid bias from un- 
certainties in the predicted absolute flux, the secondary 
calculation is normalized to the measurements at the 
spectral peak (2 GeV) as in Fig. [S] Comparing the mod- 
els shown in Fig. [3] to the measurements, only 1 and 5 
give a significant, and almost identical, excess. We use 
curve 1, a slightly better fit to the measured spectrum, to 
calculate TZ. Using the Maki et al. PBH model [5| with 
force- field modulation gives TZ = 5.0140 ^ 10~'*pc~^yr~^, 
as shown in Fig. HI This excludes by more than 9 
sigma the slight possibility of primary p's suggested by 
7^ = 4.2±1| X 10-3pc-3yr-i from BESS95-f 97 data with 
the same models and modulation. We also find a 90% 
confidence level upper limit of 7^ '--^ 1.2 x 10~3pc~3yr~^. 
This is almost insensitive to modulation (500 MV: TZ 
pc^'^yr" 

Vc""yr 

The affects of charge-sign dependent modulation in 
the A>0 and A<0 solar magnetic field polarities on sec- 
ondary and possible PBH primary p fluxes differ consid- 
erably because of their spectral shapes. Curves 2 and 6 in 
Fig. [2] indicate that the differences in the solar- minimum 



600 MV: 7^ = 1.2 X lO^^pc^^yr" 



700 MV: 7^ = 1.3 X 10" V^^yr"^) 



secondary fluxes are small. However, the predicted LIS 
PBH p spectrum peaks near the lower end of the BESS- 
Polar II energy range, and solar polarity would strongly 
affect the contribution of primaries to the measured low- 
energy flux. The primary p flux should be suppressed 
for A>0 and higher for A<0. Thus, solar polarity can- 
not explain the excess reported by BESS95+97 or the 
negligible excess in the BESS-Polar II results reported 
here. Within statistics, the BESS-Polar II data show no 
evidence of primary p 's from PBH evaporation. 
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